A problem with the escort distribution representation 
of nonextensive statistical mechanics 
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> It is pointed out that the Tsallis entropy functional represented in terms of the escort 

X 

distribution is not concave if the entropic index q is less than unity. It is emphasized 
that the escort distribution is a secondary object calculated from the basic original 
distribution. 
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In recent years, there has been growing interest in a nonextensive generalization of 
Boltzmann-Gibbs statistical mechanics. This formalism termed nonextensive statistical 
mechanics was initiated by Tsallis [1] by introducing a generalized entropy and was 
applied to various physical systems with long-range interactions, long-range memories, 
(multi)fractal structure and spatio-temporal complexity 1 . Among others, its recent 
impressive success in describing the phenomena of fully developed turbulence should 
be noticed [2,3]. 

Nonextensive statistical mechanics has a specific mathematical structure. It makes 
use of the concept of the escort distribution [4] to define expectation values of physical 
quantities. Since the escort distribution plays an essential role, some manipulations have 
been made to represent the theory completely in terms of it in the literature (which we 
do not specify here). The purpose of this short note is to call an attention to a problem 
regarding the concavity property of the entropy functional in the escort distribution 
representation. The escort distribution cannot be taken as a basic variable in maximum 
entropy calculus and is a secondary object calculated from the basic original 
distribution. 

To make the point clear, it seems appropriate to first summarize the standard 
formalism of nonextensive statistical mechanics. This theory is based on the Tsallis 



1 A comprehensive list of references is currently available at 
http://tsallis.cat.cbpf.br/biblio.htm 
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entropy 
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where {p ; } !=12L w * s me norma li ze d probability distribution, is the number of 

accessible microscopic states of the system under consideration and q is a positive 
parameter termed the index of nonextensivity . In the limit q — > 1, S q \p\ converges to 

the ordinary Boltzmann-Shannon entropy S[p] = _ E,-i P> ^ n ^ ' Difference of q from 
unity indicates the degree of nonextensivity. An important point is that S q \p\ is a 
concave functional for all positive values of q. To develop generalized canonical 
ensemble theory, we maximize the Tsallis entropy under the constraints on the 
normalization of the distribution and the expectation value of the system Hamiltonian 
H . Regarding the latter constraint, the authors of Ref. [5] introduced the concept of the 
normalized g -expectation value. Let e j be the value of the Hamiltonian in the ith 
state of the system. Then, the generalized internal energy U is given by 

w 

Z P i 9)£ i = U «> ( 2 ) 

i=l 

where P. q) is the escort distribution [4] defined as follows: 

P^ = ^—. (3) 
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The functional to be maximized is given by 



f w \ / w ^ 

*,[p;a,jB] = S,[p]-a £p,-l -jB £ /> (9) - 1/ ? 

Vi=l / Vi=l 



(4) 



where a and /? are the Lagrange multipliers. Using variational principle, this 
functional is found to be maximized by the following distribution: 
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where 
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(6) 



Z,(i8) = j;[l-(l-«)(i8/c)( e< -^)] 
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From the normalization condition on p t , follows the identical relation 
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Using eqs. (5), (6) and (8), we find another expression for Z (0) : 
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Z,(i8) = j;[l-(l-«)G8/c)(e < -^)] 
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In the limit q — > 1, the distribution in eq. (5) converges to the ordinary Boltzmann- 

Gibbs distribution p i = exp(-/?£,)/ Z(f3) with Z(f3) = £^ exp(-/?£ ( ). 

Now, let us consider the escort distribution representation. From eq. (3) and the 
normalization condition of p i , we have 



£ (p,Y = 
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Therefore, the Tsallis entropy is rewritten as 
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Now, in the escort distribtuion representation, the following functional is considered: 



f w \ ( w 
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V,[P ( « ) ;a,i8] = S 9 [P ( « ) ]-a £^-1 -j8 £^£, -[7, 
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The variation with respect to the escort distribution gives rise to 
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where Z q (f5) is the same as that given in eq. (7) [or eq. (9)]. This is precisely equal to 
the escort distribution associated with p i in eq. (5). Therefore, apparently, the escort 



distribution representation yields the same result as the standard formalism if eqs. (3) 
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and (9) are used. 

However, there are at least two points to be noticed here. First of all, the 
normalization condition on the escort distribution cannot be a constraint if the definition 

in eq. (3) is used. That is, ^ P. q) = 1 becomes an identity in this case. Therefore, if 

^J-i ^ <?) = 1 is taken as a constraint, then the relation in eq. (3) has to be discarded and 

accordingly P. q) has to be regarded as a fundamental variable. 

Secondly and more crucially, the Tsallis entropy in the representation in eq. (11) is 
seen to have a problem with the concavity condition, in general. In Fig. 1, we present 
plots of the second-order derivative, f(r; q) = d 2 S q \P {q) } j dr 2 , with respect to (r, q) 
in the case of the simplest two-state system: P^ q) = r, P 2 (9) = l-r (0<r<l). As 
clearly seen, 5 ? [P (?) ] is not concave, i.e. f(r; q) changes its sign, if < q < 1, which 
is an important regime in nonextensive statistical mechanics. Thus the theory exhibits 
instability in this regime. In this context, it may also be worth pointing out the fact that 
nonextensive statistical mechanics with < q < 1 can give rise to the negative specific 
heat [6]. 

In conclusion, we have pointed out that that the escort distribution representation in 
nonextensive statistical mechanics has a problem regarding the concavity property of 
the entropy functional and the escort distribution cannot be taken as a basic variable in 
maximum entropy calculus. The escort distribution in nonextensive statistical 
mechanics is a secondary object, which should be constructed from the basic original 

distribution {a}, =1 , 2 , l , w - 
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Fig.l Plots of f(r; q) = d 2 S q [P (q) \ / dr 2 for the two-state system, 
P^ q) = r, P^ q) = \-r (0<r<l), with respect to (r, q). 
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